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Abstract 

We consider the recently discovered Dirac semimetals with two Dirac points ±K. In the presence of elastic 
deformations each fermion propagates in a curved space, whose metric is defined by the expansion of the 
effective Hamiltonian near the Dirac point. Besides, there is the emergent electromagnetic field that is defined 
by the shift of the Dirac point. We consider the case, when the deformations are caused by the dislocations. 
The dislocation carries singular torsion and the quantized flux of emergent magnetic field. Both torsion 
singularity and emergent magnetic flux may be observed in the scattering of quasiparticles on the dislocation 
due to Stodolsky and Aharonov - Bohm effects. We discuss quantum anomalies in the quasiparticle currents in 
the presence of emergent gauge and gravitational fields and the external electromagnetic field. In particular, 
it is demonstrated, that in the presence of external electric field the quasiparticles/holes are pumped from 
vacuum along the dislocation. The appeared chiral imbalance along the dislocation drives the analogue of 
chiral magnetic effect, that is the appearance of electric current along the dislocation. 


1. Introduction 


Condensed matter systems with emergent relativistic invariance may serve as an arena for the experi¬ 
mental observation of various effects typical for the high - energy physics. In particular, some of such effects 
were predicted and checked expreimentally in 2 + 1 D graphene, where the relativistic 2 + 1 D fermions 
appear in the vicinity of the Fermi points (see the detailed review in 0 )- The emergent 3 + 1 D relativistic 
fermions appear in 3 He-A superfluid. (For the description of various emergent relativistic effects in li He 
see monograph i)- Recently, the 3 + 1 D relativistic fermions were discovered experimentally in the novel 
material - the Dirac semimetal [s, 2, 0]. For the description of the microscopic theory of Dirac semimetals 
(that may be approximated by the low energy effective field theory with emergent relativistic invariance) see, 
for example, 0- 0 , §!, 0 and references therein. The description of certain relativistic effects that may exist 
in Dirac and Weyl semimetahu was given already before the experimental discovery of Dirac semimetals 
Tol 11,12, 10 3 IS 16, 17, 3 IS IS 0H, US- There are indications that the experimental realization of 


the Weyl semimetals will also be discovered in future (see, for example, 23]). 

In this paper we proceed our previous investigation of condensed matter systems with emergent relativistic 
!6|. In the mentioned papers we concentrated on the 2 + 1 D graphene. In particular, 


m 


13, [25 


invariance 

24, 3 the emergent gravity in graphene was discussed (see also 


27 


where the similar results were 
obtained independently for the case when there are no off - plane displacements of the graphene surface). In 
flit the effects of emergent gravity and emergent gauge fields were investigated in the presence of dislocations 
(see also [3]) ■ Here we extend our consideration to the 3 + 1 D systems that simulate truly relativistic high 
energy theories. We concentrate on the Dirac semimetals. Those materials may be considered as elastic 

S3, E3, m, 


media 



and as topological matter with Weyl fermions 0 0 37. ISsf. Both Cd 3 As 2 


and NaaBi discussed in 03 0 ] have two Dirac points ±K. At each Dirac point there is the pair of Weyl 


1 Weyl semimetals differ from Dirac semimetals by the difference in the positions of the Fermi points for the left - handed 
and the right - handed fermions. In Dirac semimetals the left - handed and the right - handed Weyls fermions always come in 
pairs located at the same Dirac point. 
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fermions (the right - handed and the left - handed Weyl fermions). The chiral anomaly in Weyl and Dirac 
semimetals was discussed in a number of papers (see [39j, [ioj, [U 421 and references therein). 

Since the Dirac semimetal simulates truly relativistic high - energy theories, the discussion of anomalies 
in the latter theories may be applied to it (with some reservations, though). In this respect it is worth 
mentioning that there is a certain contradiction in the results published to date. Namely, it was reported 
that in the presence of torsion chiral anomaly includes the topological Nieh - Yan term. In particular, this 
has been demonstrated using Fujikawa method (43[. In |44j| . however, it was argued (this was based on the 
perturbative calculations with Pauli - Villars regularization) that there is no Nieh - Yan term term in the 
chiral anomaly. In [40j the alternative derivation of the Nieh - Yan term in the anomaly was given that is 
based on the consideration of massless fermions of opposite chirality as living on the boundary of the 4+1 
D topological isolator. In the present paper we do not discuss the situation, when the Nieh - Yan term may 
be present in chiral anomaly. Instead we restrict ourselves by the case of constant in time vielbein with 
vanishing components e l 0 (and vanishing spin connection), in which the Nieh - Yan term vanishes. 

We derive expressions for anomalies in quasiparticle currents using the explicit solution of Pauli equation 
in the presence of nontrivial vielbein that does not depend on time. The spectral flow in the presence 
of gauge fields and the torsion magnetic field results in the appearance of the anomalies in quasiparticle 
currents. Notice, that this method was also applied in [40] . However, the emergent magnetic field was 
not taken into account (in case, when the torsion magnetic field is present). At the same time it always 
accompanies the emergent gravity in Dirac semimetal. As a result the pattern of the spectral flow and the 
pumping of quasiparticles from vacuum that we discuss in the present paper differs somehow from that of 

The emergent gauge and gravitational fields experienced by fermionic quasiparticles in Dirac semimetal 
in the presence of elastic deformations are expressed linearly in those deformations. The coefficients entering 
these expressions generalize the so - called Gruneisen parameter of graphene 3, 251. The elastic deformations 
in the presence of dislocations may be calculated using elasticity theory |45^ . It is worth mentioning, that 
in the presence of dislocations the emergent vielbein and emergent gauge fields are nontrivial even when the 
contributions of the material parameters is neglected (see also [46j). In the present paper we assume, that the 
material parameters may indeed be neglected in a certain approximation, and consider the purely geometrical 
contributions to emergent vielbein and emergent gauge field. The spin connection in this approximation 
vanishes. For the early discussion of the concept of emergent gravity in the system of interacting relativistic 
fermions see, for example, [47l ] and references therein. 

According to our results the dislocation in Dirac semimetal as well as in graphene carries both emergent 
magnetic flux and torsion singularity. The torsion singularity may be observed in the scattering of the 
quasipartcles on the dislocation and results in Stodolsky effect. The history of this effect originated from 
the suggestion of Colella and Overhauser [3] to measure the beam phase shift in the real gravitational field 
(that is expressed through both gravitational constant and Plank constant). This effect was also discussed 
in a number of papers (see, for example, [ 49 J). Stodolsky was the first who derived [i5C| an expression for the 
phase shift in the case of the gravitational field of general (non - Newtonian) form. The observation of this 
effect in graphene for the scattering of quasi-particle on the torsion singularity carried by the dislocation was 
proposed in 129]. 

Finally, we calculate various effects of quantum anomalies in the quasiparticle currents in Dirac semimetals 
and apply our general expressions to the case, when the dislocations are present. It is worth mentioning 
that the anomalous production of the quasiparticles does not mean that they appear from nowhere. The 
microscopic currents are conserved, and the appeared charge refers to the charge of the free quasiparticles 
only. Suppose, that the left - handed electron at the Dirac point K appears at the dislocation (in the 
presence of electric field). We may consider this process as a jump of the left - handed electron from one 
zone (electron sea in vacuum) to another (conductance zone, the zone of quasiparticles). The total number 
of electrons of each chirality remains the same, but the number of the left - handed cells in vacuum at the 
dislocation is decreased. Therefore, the electron that occupied a cell in vacuum before the appearance of 
electric field, becomes free when the electric field is turned on. 

Pumping of quasiparticles from vacuum results in the appearance of chiral imbalance along the dislocation. 
This drives the analogue of the chiral magnetic effect that is the a ppe arance of electric current along the 
dislocation. Notice, that the similar phenomenon was discussed in 15 ll |. where the external magnetic field 
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plays the same role as the emergent magnetic field of our case. 

It is worth mentioning, that in the present paper we neglect completely Coulomb interactions between 
the quasiparticles. These interactions are strong due to the large effective coupling constant [55| and break 
the emergent relativistic invariance. However, the investigation of graphene indicates, that such interactions 
may be neglected for the consideration of certain properties of the system [lj. We expect, that the Aharonov 
- Bohm and Stodolsky effects in the scattering of the quasiparticles on the dislocations as well as the chiral 
anomaly driven by the emergent magnetic field incident at the dislocation are such effects that survive when 
the Coulomb interaction is turned on. 

The paper is organized as follows. In Section [2] we briefly discuss the general theory that describes 
reduction of multi - fermion systems to the low energy effective models with Weyl fermions. In Section [3] we 
discuss the appearance of emergent gauge fields and emergent gravity with torsion experienced by fermionic 
quasiparticles in the Dirac semimetals in the presence of elastic deformations. In Section [4] we calculate 
the emergent magnetic flux and the emergent torsion carried by the dislocation (assuming that the material 
parameters may be neglected). We also discuss how the emergent magnetic field and the torsion singularity 
may be observed via the Stodolsky and Aharonov - Bohm effect. In Section[5]we consider quantum anomalies 
in quasiparticle currents in Dirac semimetal, and apply the obtained general expressions for the calculation 
of the anomalous production of quasiparticles at the dislocations in the presence of external electric field. 
Also we discuss chiral magnetic effect caused by emergent magnetic field. In Section [G] we end with the 
conclusions. 


2. Reduction of multi - fermion systems to the models with Weyl spinors. 


2.1. Microscopic theory 

In this section we briefly discuss reduction of a multi - fermion system of general type to the model with 
Weyl fermions. This subject has been discussed in a number of papers (see [2| and references therein). The 
general theorem that describes this reduction has been proposed in [37]] ■ It may be applied to the wide class 
of condensed matter systems with Fermi points. The Dirac semimetal is the particular example of such a 
system. In this section we follow the narration of 
Let us denote the original n ■ 
the form: 


component spinors by tjj. The partition function of microscopic theory has 


Z = J DipDijjD^expyiR[^] + i 




( 1 ) 


The collection of fields that provide interaction between the fermions is denoted by $. I? is a certain function 
of these fields. In mean field approximation the values of $ are set to their ’’mean” values $o and we omit 
integral over $, and denote H(V) = H(§ 0 , V). )T) X may be understood as the integral over d 3 x in continuous 
systems and may be considered as the sum in discrete systems. 

By V let us denote the three commuting Hermitian operators, whose eigenvalues parametrize the branches 
of the spectrum of H . Those parameters do not necessarily coincide with the ordinary three - momentum. 
In the basis of the corresponding eigenvectors the fermion field is denoted by ift-p. By construction the 
coefficients in the expansion of H in powers of V do not depend on coordinates. Different branches of 
spectrum of H repel each other. This means that any small perturbation pushes apart the two crossed 
branches. That’s why only the minimal number of branches of spectrum may cross each other. This minimal 
number is fixed by momentum space topology Q- 

Let us suppose, that at V = V 1°' the n re duced branches of H cross each other. Then there exists the 
Hermitian matrix H( V) such that H(V) = (V)H(V)fl(V) is diagonal. Its first ?r re duced x ?r re duced block 
-^reduced corresponds to the crossed branches and to the fermion components 'Fp = = 1,..., n re duced- 

All eigenvalues of H le duced('U) coincide at V = V The remaining block H ma ss i ve corresponds to the 
’’massive” branches and to the fermion components Q J V = -0p +nreduced , j = 1 — n re duced- The functional 

integral can be represented as the product of the functional integral over ’’massive” modes and the integral 
over n re duced reduced fermion components. 

Using transformation of the form ip-p —>• e~ tEot ipj> , H{V) —> H{V) — Eq (that means an adjustment of 
chemical potential) we can always make the only eigenvalue of iL r educed(T , *°' ) ) equal to zero. Therefore, V ^ 
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may be considered as the position of Fermi point. At the same time variable V is to be treated as generalized 
momentum. As a result the partition function has the form 


Z = 


D'i’D'i’DQDQexTpli 


dt^2\^ V {t)(idt - ^reduced + 0 V (t)(id t - ^massive )&v(t) ) 
V 


( 2 ) 


2.2. Reduced low energy model 

The eigenfunctions of iA-educed do not depend on time at the Fermi point. Therefore at low energies 
the integral over T-p dominates. The other components contribute the physical quantities with the fast 
oscillating factors. As a result they may be neglected in the description of the low energy dynamics. As it 
was mentioned above, the branches of spectrum repel each other. Therefore, the branches crossing occurs 
only if it is protected by momentum space topology. The minimal number of fermion components that 
admits nontrivial topology is two Q (for the discussion of topological invariants see the next subsection). 
Then the reduced Flamiltonian has the form 

^reduced = M%(V)a k + M(P) (3) 

Here M^(V) and M{V) are the real - valued functions of the generalized momentum V. When operator 
V is substituted instead of real - valued vector we arrive at the Hermitian operators M k (V) and M(V). 
Expression of Eq. d3j) is to be understood as the Hamiltonian in momentum space, i.e. in the basis of the 
reduced fermions (the eigenvectors of the generalized momentum are the elements of this basis). Our next 
purpose is to represent this Hamiltonian in the long - wave approximation in the basis of the eigenvectors of 
coordinates (that is we need to rewrite Eq. ([3]) in coordinate representation). The explicit solution of this 
problem is not so evident in the case, when V are Hermitian operators of general type. The eigenvectors of 
these operators that correspond to the crossed branches are expressed as complicated linear combinations of 
the basis vectors of coordinate representation. 

According to our notations space of 'Fp is spanned on the eigenvectors | V,i), i = 1,2 (while space of 
massive components is spanned on the eigenvectors \P,i), i = 3, ...,n), at the same time the eigenvectors 
of coordinates are |x, i), i = 1 , ...,n and we may write symbolically (x, i\i/j) = ■?/»*( x ) and (V,i\ip) = ( we 
should not forget, however, that ij) is Grassmann variable). Our aim is to construct the linear combination 
of | V,i), i = 1,2 in such a way, that it is the (approximate) eigenvector of the coordinate operator x. We 
accept the natural assumption, that the definition of vectors |x, i)’ exists such that: 

1. (y, 7 lx, i)' —> 0 at lx — y| 1/A for some scale A; 

2- (y, j| x , i)' ->• 1 at |x - y| <C 1/A; 

3. |x,*)' = X]-p;j=i ,2 Qij('P’ x )\'P’j) f° r bi = 1,2 with some functions qij(V,x). 

In the other words, we may compose the wave packets of linear size ~ 1/A (localized near any given point 
in space) out of the eigenvectors of V corresponding to i = 1,2 (that is out of the states with definite values 
of energy corresponding to the two crossed branches). The uncertainty principle tells, that this may not 
be done only if the states with definite values of energy (corresponding to the two crossed branches) have 
distinct discrete values of at least one component of the real momentum p (its eigenvectors have the form 
~ e px ). We assume that our system has an infinite size and do not consider here such a marginal case. 

Thus for the energies below a certain value A we may define 

= <x,i|V) = ^2 * = 1,2 (4) 

r-,j= 1,2 

In this approximation we replace summation over coordinates by an integral, and write the reduced Hamil¬ 
tonian as 

^reduced ~ Wfc (x, p)<J fc + m(x, p), (5) 

where p = — iV is the ordinary momentum in coordinate space. Functions to^,to are Hermitian operators 
that depend both on momentum and on coordinates. Eq. © should be understood as 

J2^r(t)(Mk(P)a k + M(P))V-p(t) « f d 3 x§(x, t)(m%(x, p)a k + m(x, p))^(x,t) ( 6 ) 

-p J 
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In the particular case, when , m do not depend on coordinates the generalized momentum coincides with 
p = — iV, while M k (V) = m k (fP), and M(V) = m{V). 


2.3. The role of momentum space topology 
The partition function is reduced to 


Z = 


DU/DU/expf* / d‘ l x'S/(x)(idt — m k (x, p)<r fc — m{x 1 p))\E'(a; 


(7) 


We assume, that the dependence of m k ,m on coordinates is weak and therefore in the following we may 
formally deal with functions m as independent of coordinates. This allows to introduce the notion of the 
floating Fermi point K(x) as the value of p at which m k (x, p)a k vanishes. The value of K also contains slight 
dependence on coordinates. When the dependence of Hamiltonian on x is absent the value of K coincides 
with the position of the true Fermi point, and m(K) = 0 as well because in this case variable p coincides 
with V. 

The energy spectrum is E{V) = ±\M L (V)\ + M{V). If all three components of V are continuous, the 
topological invariant that protects branches crossing has the form 


n = f dS l M L 
8 t r J a 


f dM L dM L \ 

\dv~ x J 


M L 


M l 

\M L \ 


( 8 ) 


where cr is the S 2 surface around the point . For N = 1 the branches crossing cannot disappear as a 
result of the continuous deformation of the system because the zero of A4jt(P) cannot disappear while the 
second term M(V) in Eq. ([3]) gives the same correction to both eigenvalues of ^reduced- In this case the 
expansion of the functions M k {V) and Af(V) near the Fermi point has the form Q: 

m? op) » 

M{V) » Eo + Ji{V j -V? ) ) (9) 


Here J-j are the real - valued coefficients of expansion. The quantity Eo may always be eliminated by the 
adjustment of chemical potential. The case N = 2 would correspond to the expansion, in which the term of 
the first order in (Vj — V^) is absent in the expansion. The configuration of Aif(V) with a certain value 
of N cannot be transformed continuously to the configuration with different value of N. 

In the situation, when one or more components of V may be quantized the consideration is more involved. 
It is supposed, that our system of general type is obtained as a continuous deformation of a system with mjt 
and m independent of coordinates (in which case those functions coincide with M.^(V) and A4(V)) with the 
topological invariant of Eq. © equal to N = 1. In this case we are able to expand operators mf and m in 
powers of p near the floating Fermi point : 


mf(x, p) 
m(x, p) 


l - K,-(x)) + (pj - K Ax))fl (x) 


2 i 


1 


A 'o(x) + - fl (x)(pj - Kj(x)) + (p j - Kj(x))/q (x) 


( 10 ) 


Here we take into account that mf and m are Hermitian, which results in the symmetric form of the product. 
By f[ (x) and Kq(x) we denote the coefficients of expansion. If their dependence on x is neglected we would 
arrive at f- « J 7 ] and K 0 ss 0. 


2-4- Emergent Weyl spinors 
Now Eq. © has the form: 


Z= -D’FD’Fexpfi- 


dfx 


^{x){id t - K 0 (x) - f 3 k {x){f>j - K,)6- fe )^(x) + (h.c.) 


( 11 ) 
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Here ( h.c .) means hermitian conjugated expression. According to our definition the hermitian conjugation 
applied to dr gives T. In the partition function of Eq. ED, the sum is over k = 0,1, 2, 3, and j = 1,2,3 
while cr° = 1. 

The coordinate dependent matrix /®(x) coincides with the effective vielbein e„(x) up to the factor 
e(x) = det 1 / 3 ^/*(x)^ that defines invariant integration measure over coordinates: 

/i(x)=e(x)e*(x) (12) 


We require e° = 0 for a = 1, 2,3, and e x e[j = 1. Here e _1 = e[j x det 3 X 3 e\ = e q is equal to the determinant 
of the vierbein e®. We arrive at the partition function of the low energy effective model with emergent 
relativistic invariance 


Z = 




(13) 


where 


S 


1 

2 


d 4 x e 'f>(x)e J a d a i£)j'fy(x) + (. h.c .)^. 


(14) 


Here the sum is over a,j = 0,1, 2, 3 while a° = 1, and D is the covariant derivative that includes the [7(1) 
gauge field K fl = (K 0 , K). 

If we want to take into account interaction between the fermions due to the fluctuations of the fields $ 
in Eq. ©, then we should start from Eq. © with the fixed value of $ that is not equal to its mean value 
$ 0 . At each value of <I> we proceed through the above steps and come to Eqs. (fl3l) . (ThU) with the values 
of e J a and K fl dependent on $. This dependence may be rather complicated. Nevertheless, in principle, the 
expansion in powers of $ — $o with phenomenological coefficients may be used. This will bring us to the 
theory of relativistic fermions interacting with the collection of bosonic fields. 

One can see, that in the considered long wave approximation the emergent spin connection does not 
arise. Actually, we may add the spin connection for completeness. But for the Hermitian Hamiltonian it will 
effectively result in the interaction of the Weyl fermion with a new [7(1) gauge field and to the correction to 
chemical potential. Thus, those contributions may be absorbed by the field K M and by the overall chemical 
potential. 


3. Emergent Weyl fermions in Dirac semimetal in the presence of elastic deformation. 

3.1. Dirac semimetals as elastic media 

First of all let us consider the Dirac semimetal as elastic media. The deformations in 3 D elastic media 
are described by the change in the coordinates of atoms X‘(x) = u l (x) + x l , where u l (x) is the displacement 
field, which describes the displacement of atoms from their equilibrium positions. We assume that the elastic 
deformations are small diU k -C 1. 

The displacement field is expressed in terms of the elastic deformations and rotations: 


djUi(x) = tij(x) + u)ij(x) (15) 

C-ij ~ 2 + 71/ a,j ), c Jij = — (0j '/.// diUj) (Id) 


Tensor uiij describes local rotations, while tensor of elastic deformations is denoted by e-ij. Throughout the 
text we assume that elastic deformations do not depend on time. 

Elasticity theory operates with metric 


dX l dX l 
dx l dx k 


dik + 2 Cik 


(17) 


The naive supposition is that the fermions propagate in Riemannian space corresponding to this metric. 
However, we will see below, that as well as in graphene 24, [25|, the fermionic quasiparticles in Dirac 
semimetal in the presence of elastic deformations propagate in Riemann - Cartan space with vanishing 
curvature (Weitzenbock space), which is determined by a certain vierbein. In leading approximation this 
vierbein is expressed linearly through . 
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3.2. Foui ’ Weyl fermions in Dirac semimetal. 

It follows from the microscopic theory [2, 0, 0, 0] that in the considered Dirac semimetals (without strain) 
there are four Weyl fermions aaa (the pair left - handed fermion right - handed fermion at K^ 0 -* and the 
pair left - handed fermion - right - handed fermion at — K7 0 )). Thus in the absence of elastic deformations 
we have two pairs of coinciding Weyl points ±K7 0 ) and the low energy effective action of the form 

S = i ^2 Y, [ d3 xdteW 9 St ±e%K?(idv + KW ± )* Bt± + (h.c.l a° = 1 (18) 

s=L,R ± •' 


where <7p = a K while erf = <j k , and o a = —o a for a = 1,2,3, while o° = er 0 = 1. dq 'b are the independent 
fermion Grassmann variables that describe quasi - particles living near the Weyl points ±K^°fx). Here 
Kl 0)± = ±K™> for k = 1,2,3 and A^ 0)± = 0 while 


= Vp ; 


o(0)i — 

-'a J ai 


AO)i 


= 0 , 


p(0)0 _ 

'-'n 


VF 


_A° 

°0> 


= 1,2,3 


(19) 


where vpfa has the meaning of anisotropic Fermi velocity. The considered Dirac semimetals 0, Q 0/ are 
anisotropic, that results in the anisotropic Fermi velocity corresponding to 3x3 matrix f = diag(zz _1//3 , z/ -1 / 3 , z/ 2 / 3 ) 
with z/^1. 

We do not exclude, that in the presence of elastic deformations the positions of the two Dirac points 
may split into four (floating) Weyl points depending on the position in space. We denote them by K L,± 
and K fl,± . The definition of K depends on the choice of coordinate system, i.e. on the parametrization 
of the semimetal. We may choose the parametrization, in which the coordinates of atoms in the crystal 
are the same as in the unperturbed lattice. In this reference frame we have K L,± ss -f- A(“) l,± and 

K fli± « ±K<°) + A(“) r ' ± , where A^’* and Aare to be interpreted as the emergent 17(1) gauge 
fields. We refer to the given parametrization as to the accompanying reference frame. 

If the crystal is deformed only slightly, the other parametrization is preferred, which is called typically the 
laboratory reference frame. In this parametrization the coordinates of atoms are their real 3D coordinates. 
Then the transformation between the two parametrizations is given by: 


X k (x) = u k {x)+x k , k= 1,2,3 


( 20 ) 


Here the laboratory reference frame coordinates are denoted by X k while the coordinates of the accompanying 
reference frame are denoted by x k . 

We neglect possible interactions between the fermions located at different Weyl points. Therefore, the 
general theory described in the previous section may be applied to each of the mentioned Weyl fermions. As 
a result the effective low energy action has the form of the sum over the four Weyl points: 

S = i ^2 Y, ( d 3 xdt e 1 Fs,±e^ =l= ’ At of : (id^ + \1G,± + ( h.c .), er° = 1 (21) 

L s=L,R ± ’ 


d', d' are the independent fermion Grassmann variables that describe quasi - particles living near the (floating) 
Weyl points K s,± (x). The effective 3 + 1 D vielbein may be different for the four Weyl fermions. 

Therefore, we add the superscript s,± . Here K^ for k = 1,2,3, while may be nonzero. 


3.3. Expressions for emergent vielbein and emergent gauge field. 

For small elastic deformations the values of emergent gauge field and the values of the vielbein should 
be expressed through the derivatives of the displacement vector u k . In the decomposition of Eq. (USD the 
contribution of ujij corresponds to rotations that cannot cause strain and therefore does not contribute to 
the emergent vielbein and to emergent gauge field in accompanying reference frame. In laboratory reference 
frame ujjj may appear, however. The explanation of this difference follows the analogy with graphene 24, 25( | . 
where it results from the consideration of microscopic theory. Explicitly, in Dirac semimetals as well as in 
graphene the microscopic model operates with the fermion operators ip(x) defined on the sites of the lattice 
(see 0,0, ii and references therein). These operators create and annihilate electrons on certain orbitals 
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that belong to the atoms located at the lattice sites. The microscopic Hamiltonian depends on real distances 
between the atoms through its phenomenological parameters. Small elastic deformations result in small 
variations of these parameters. 

The accompanying reference frame is distinguished because in this reference frame the coordinates of 
atoms are the same as in the unperturbed lattice. This is this reference frame, which appears naturally in 
the low energy effective continuous model because operators i/>(x) are attached to the sites of the lattice. 
Therefore, this is this reference frame, in which the parameters of effective action Eq. HID e v a and K M are 
equal to their unperturbed values plus corrections that are proportional to the derivative of the displacement 
vector u l . As it was mentioned above, tensor ujij in the decomposition of Eq. m corresponds to rotations. 
Rotations of crystal cannot change the distances between the atoms. Therefore, uj i:) does not enter the 
expressions for e v a and K^ in this reference frame. For the mentioned above reasons this consideration refers 
to the accompanying reference frame only and does not work, for example, in laboratory reference frame, in 
which the coordinates of atoms are given by their real 3 D coordinates. In laboratory reference frame the 
rotation of crystal results in rotation of Fermi point K (while in accompanying reference frame rotations 
cannot change K). In the same way in laboratory reference frame rotations of crystal result in rotation of 
vector e l a for each a = 0,1, 2,3. It will be seen further, that in Laboratory reference frame tensor enters 
expressions for the emergent vielbein and emergent gauge field due to the transformation between the two 
reference frames, which contains both Wy and e,- L j. 

Thus, the most general leading linear dependence of gauge field on e,j in accompanying reference frame 
has the form 


A (.a)L,R,± 


1 oL,R,± A {a)L,R,± 

Pijk SIq 



i,j,k = 1,2,3 


( 22 ) 


where = ±K^ + a[“) L ’' R ’ ± for k = 1,2,3 and Kq' r,± = A^ L,R,± , while is the tensor of 

dimensionless material parameters (similar to the Gruneisen parameter of graphene), a is the lattice spacing. 
In laboratory reference frame the values of K are given by: 


K 


dx h 

dX‘ 


( 


±K 


(o) 


(o)± 


) « (s* - diU k ) ( 


±K 


( 0 ) 


, (fo± 


(±K-°^ V;(u • K^)) + A- a)± 


(23) 


Here we omit the superscripts L,R. The field A* 0 ) is given by Eq. (l22l) . In the following we denote the 
emergent gauge field in laboratory reference frame by 


A f » tV ! ;(u-KW)) + A 1 (q)± 


(24) 


The first term in the right hand side of this expression contains the geometric contribution, which comes 
from the coordinate transformation of the original position of the Dirac point in the non-deformed 
lattice. 

In the same way we relate parameters /* = ee l a of the fermionic spectrum with the elastic deformation. 
As well as for K the definition of /£ depends on the choice of coordinate system, i.e. on the parametrization 
of the semimetal. In the accompanying reference frame we have: 

fa ~ V F (Jl - lajfeGfc) , /o ~ -VFlbjktjk, fa = C <b i, j, k = 1,2,3 (25) 

where at each Weyl point the tensor of new dimensionless material parameters 7’ is defined (we omit here 
the superscripts L , R, ±). 

Applying coordinate transformation to Eq. (12511 we come to the following expression for /* in the 
laboratory reference frame: 

fa = det" 1 (UK (fa - 7 Zj&l) U * v F (/* (1 - d k u k ) + tdkvf - 7 Zjftl) (26) 

Then we come to 

fa ~ V F (f a ( 1 - e kk ) + f k e ia + f k UJia - 7 - diUj ) (27) 
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The terms /*( 1 — dkUk) + f k dkU 1 in the right hand side of Eg. ESI) are of the geometric origin: they come 
from the coordinate transformation from the original non-disturbed to the strained semimetal. 

We come to the conclusion, that when the material parameters /3ijk and 7 * are neglected, in laboratory 
reference frame we arrive at the same expression for the vielbein for all four Weyl fermions that is given by 
the pure geometric contribution. In this approximation (which will be assumed in the following) the spin 
connection vanishes while the emergent gauge field is again the same for all four Weyl fermions, is given by 
the first term of Eq. E4l) and has opposite charges for opposite Dirac points iK^ 0 ). 

The contributions to various physical quantities of the material parameters may be neglected if those 
dimensionless parameters are small. In the following we assume that for the considered Dirac semimetals 
the values of /3ijk and 7 * . fc are indeed small. This means that elastic deformations with e,^ ~ 1 result in the 
contributions to the quasiparticle Hamiltonian that are small compared to the unperturbed Hamiltonian. 
More explicitly, first, the elastic deformations should result in the contribution to the emergent gauge field 
that is small compared to its unperturbed value K^. This is a very natural assumption because at the 
same time it is required in order to apply low energy effective field theory. (It is supposed, that this theory 
may be applied if values of momenta including its shift due to the emergent gauge field are small compared 
to K^°i.) Second, elastic deformations should give contributions to the anisotropic Fermi velocity that are 
small compared to its unperturbed values. Qualitatively this is also natural, because this means that elastic 
deformations are not able to change considerably speed of various excitations existing in the semimetal. 
Nevertheless, the mentioned assumption is to be checked experimentally for the particular materials. 

4. Emergent geometry in the presence of dislocations 

f.l. Accompanying reference frame 

Effective geometry experienced by fermions acquire singular contributions to torsion concentrated at the 
origin 1° of the dislocation (that is the closed line or a line connecting the boundaries). At the origin l 0 of 
the dislocation the cells of the crystal may be substituted by irregular cells. Then the extra sequence of 
regular cells is added along the surface J (line l 0 is its boundary). For the low energy effective theory this 
results in cutting of the volume along the surface J that begins at 1° and goes to infinity. Then the strip of a 
finite width is added along J. The resulting volume is sewn along the cut. As a result in the accompanying 
reference frame (where the coordinates of the atoms are the same as in the unperturbed lattice) there is the 
uncertainty in the definition of the parametrization at the cut. This uncertainty gives nonzero value to the 
following integral along the contour C surrounding 1°: 

b* = J dx l (28) 


f.2. Laboratory reference frame 

In laboratory reference frame the parametrization of semimetal volume does not contain ambiguity. The 
ambiguity in the parametrization of accompanying reference frame appears through the displacement field 
defined as a function of X k : 

x k = X k - u k (X) (29) 

In accompanying reference frame the parametrization is defined modulo the step - like discontinuity at the 
cut. Therefore, although the displacement field u l (X) has a step - like discontinuity concentrated along the 
cut, its derivative djU k is continuous. The Burgers vector is 


b l = dx l = - du l 


(30) 


For 1° of the form of the infinite straight line we may choose, for example, the following representation: 


b a 

U a = + < 0 „t: 


(31) 
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where <j> is the polar angle in the plane orthogonal to the dislocatiorU while tt“ ont is continuous along the cut 
(i.e. it does not contain the discontinuity along the cut). However, ii con t may be undefined at the origin of 
the dislocation. In the following we will use the laboratory reference frame, where there is no ambiguity in 
parametrization. 


4-3. Elasticity equations 

Elastic part of the free energy may be written as |45| : 

F=— J d 2 X + 2 /ieifcCifc^, (32) 

where A and fj, are Lame coefficients. We may neglect the term quadratic in u a in e^. Therefore, the free 
energy is quadratic in u a (a = 1,2). As a result the differential equations are linear in u a . 

The variation over u a (a = 1, 2) gives: 

d k £ik + 1 fT r) diCu = 0, (33) 

1 — AG 


Here we introduced Poisson parameter a as: 1 ff 2 a = fi ■ ^q. (ESI together with the boundary conditions 
determines the values of the equilibrium displacement field u a . Let us consider the case, when 1° is infinite 
straight line. The problem becomes two - dimensional. Its solution is given in (45|. For the completeness we 
present it in Appendix A. The resulting value of u k ont (X) is given by 


t(X) 


1 - 2(7 

2 ( 1 -^) 


,, v _L, og Ai£: 

2tt B 2 R 


+ b l i il 


1 X j X k 
1 — cr 4-7T 


(34) 


This expression should be regularized both at |A'| —> 0 and at |X| -» oo. The considered effective field 
theory fails at the distances of the order of the lattice spacing a. The effects of the finite size of the sample 
are also strong. It is worth mentioning, that while the values of u k ont may be large, its derivatives are small 
for sufficiently small b because after the differentiation the expression in Eq. (HH 1 ) tends to zero at |X| —> oo. 


4-4- Emergent magnetic field carried by the dislocation 

In order to calculate the emergent magnetic field we use integral expression 

i ijk [ H l dx J A dx k = [ A k dX k (35) 

JS JdS 

For the considered solution of elasticity equations we represent the right hand side of this expression as 
follows 

[ A k dX k = ± 6 ‘K[ 0) + 0(/3) (36) 

JdS 

The first term in this expression gives the following singular contribution to magnetic field: 

H k ng « ±b i K™ f dY k (s)5&(X - y(s)), (37) 

The unperturbed Fermi point is defined up to the transformation —>• K(°) + G, where G is the vector 
of inverse lattice. This corresponds to the change of the magnetic flux by A<I> = b • G = 2irN, where N is 
integer. Such a change of the magnetic flux is unobservable for Weyl fermions. 

According to the remark presented at the end of Section [3l we assume, that the mentioned singular con¬ 
tribution to magnetic flux dominates over the regular one proportional to parameters f3ij k - That means that 
the dimensionless generalized Gruneisen parameters (3ij k are small for all four Weyl points. This supposition 
is also justified by the analogy with graphene. We also suppose, that the dimensionless parameters 7 ® 
are small, and the corresponding contributions to torsion may be neglected compared to the singular one 
considered in the next subsection. 


2 if the third axis is directed along the dislocation, then X 1 = X® + r cos 0. X2 = A'!) + r sin 0 
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f.5. Torsion carried by the dislocation 

In the absence of dislocation the torsion tensor is defined as 


rpa 

1 jk 


fi e k \1 


(38) 


where e k is the (dimensionless) inverse 3 + 1 D vielbein related to /* according to Eq. (1121) . In the presence 
of dislocation we use integral representation: 


1 

2 


T“dX l A dX j 

L J 


[ e a k (X)dX k 
J as 


(39) 


In order to calculate torsion at A' 0 we should choose S as its small vicinity. 

In the accompanying reference frame at each Weyl point for the components of e k with k, a = 1, 2, 3 we 
have: 


e%{x) = (/*- 


0 ( 7 ), 


■fct 1/3 (/, 

r 1 


a - 7 ajk^jk) 

= diag(i/ 1 / 3 , i/ 1 / 3 , p~ 2 / 3 ) 


The values of e k in laboratory reference frame are given by 


(40) 


e a k {X) = g^e“(x) « (r a k -r?d k ul + 0( 7 )) (41) 

In principle, we should also transform derivatives entering e ka '- gfv —t However, we are able to 

substitute here jfpr by S k in the approximation linear in displacement u z ,i = 1,2,3. 

According to Eq. (l39l) torsion is related to the circulation of the inverse vielbein e k along the closed 
contour. For the case, when dkU a diu a (a = 1, 2, 3) may be neglected we have: 


J e k (X)dX k » J (-f“S fc ^(X) + 0( 7 ))dX fc = f“6 J + 0 ( 7 ) « 6“ (42) 

Recall that f k = diag(zA/ 3 , z/ 1 / 3 , v~ 2 / 3 ) and f k = diag(i/ -1 / 3 , 1 z -1 / 3 , z/ 2 / 3 ). Then b a = f k b k . We arrive at the 
following singular contribution to torsion: 


T%, sing « b a i ijk f dY k (s)S^(X-Y(s)) + 0{ 7 ) 


(43) 


^.£>. Observation of torsion singularity and emergent magnetic flux 

Here and below we use the following equivalent notations for the four vectors: X M = (A°, X 1 , A 2 , A 3 ) = 
(t, r) . Let us estimate the effect of torsion singularity on the scattering of quasiparticles on the dislocation. 
The wave function satisfies Pauli equation 

l<T?, R (ee&iD li + iD li et&)* = 0, //, K = 0,1,2,3 <t° = 1 (44) 

The covariant derivative D M contains both emergent magnetic field and emergent spin connection. Recall, 
that the vielbein is given by: 

< = k + ftd v u» + 0( 7 ), a, fi = 1,2,3 # = —(l + 0( 7 ))# 

vf 

e °ii = ^- 7 ^ 0 ^" + 0 ( 7 ), a, 11 = 1,2,3 e° =u f (1 + 0( 7 ))J° 
e = det e“ = uf(1 — dtu 1 + 0( 7 )), i = 1,2,3 (45) 


In the following we neglect all terms linear in material parameters 7 , /?. Let us represent 

k '< (Y ) dYM+i f(t 0 l) A ^ Y)dY " 


+ (r, t) = +k(r, t)e ^ * 0 ' r o> 


(46) 


ll 





where k y = (ko,k x ,k yj k z ). We omit the superscripts ± that mark the valleys for simplicity. The integral 
here is along the trajectory , given by the function Y p (s, X). It is parametrized by s and depends on 

the endpoint X p = (i,r). We assume, that is slow varying. Then vector k a plays the role of the three - 
momentum of the incoming particle because around the point (foXo) situated far from the dislocation the 
vielbein e“ « diag(qp, zA/ 3 , z/ 1 / 3 , j/ -2 / 3 ) is flat, and 


vjf ^ e -TOFfeo(t-io)+iPj(r J -rJ) 


(47) 


Here the eigenvalue of the 3 - momentum —zV is denoted by p., = —fjk a . Eq. m actually defines a 
series of different solutions that correspond to the difference in a winding of C around the dislocation. Phase 
i & a e“(y)dE Al in Eq. (PTOT) may be interpreted as follows: the translation between the two points 

X 0 = (*o,r 0 ) and X = (f, r) in Weitzenbock space is defined as R a = e^ j (Y)dY IJ '. Plane wave is 

determined by phase k a R a . Substituting Eq. (H7>1) to Eg. (HU) one obtains 



dY u (s, X) 


dY p + iV L 


)^k(r,t) = 0. 


(48) 


Here we denote T M = ^e“Vi,ee)( = 

We can always choose the trajectories Y(s,X) far from the origin of the dislocation (at the distance 
|r| » a.) For the vielbein given by Eq. (TUT) with u a given by Eq. (l34l) we have V^ee^ ~ The non - 

singular part of torsion is also ~ Therefore, is constant and satisfies aL, R ki'&k = 0. 

The infinite straight line 1° of the dislocation plays the role of a line - like hole in 3 D volume of the 
semimetal (or a 2D surface in 3 + 1 D space - time) carrying singularity of torsion Xj a 2 = b a + 0(^)6^ (X). 
The difference between the phases of the two solutions defined by paths and ended at the same point 
X M is defined by the winding number Kij of the contour C^z) = C 1 ' 1 ' 1 — around 1°. The two solutions 
differ by the phase factor: 




W 


(X) = 'icM ejfcadx^+t J c(y) Ay.dX* _ 


p iK ij b k K k +0{ t,/3,A) 
^ k c 1 


Kk=± K 


( 0 ) 


P k 


(49) 


Here /C may treated as the total momentum of the quasiparticle (that is the sum of the unperturbed Dirac 
point ±k7°) and p). Following [26[ we call the appearance of this phase in the wave function (proportional 
to the the winding number and to the momentum of particle) the Stodolsky effect. The observation of this 
effect in graphene was first proposed in [2_9fl. Actually, our consideration of the present section follows the 
approach of j2f| in its form developed in |2(| . 


5. Effects of anomaly in the presence of dislocations 


5.1. General expression for the anomaly 

The partition function of the effective low energy field theory of Dirac semimetal has the form: 


Z= j DTDTexp(*i J2 J2 [ 


s—L.R ± 


e s ’ ± 4' Sl± e^ ± ’Vf(*5 /J 


— B y ) \& s> ± + ( h.c .) 


(50) 


We introduced here the external 17(1) gauge field B, take into account that K’R = ±K^ +A* ,± and assume 
that ef k diag(l/uF, z' -1 / 3 , 1 / 2 / 3 ). 

We perform the following coordinate transformation: 

X° Vp'X 0 , X 1 ’ 2 ->■ v~ 1/3 X 1,2 , X 3 -)• v 2/3 X 3 (51) 

and the corresponding transformation of the other vectors and tensors (with the space - time indices /q v ,...) 
including K y and u p . Let us neglect all effects that result from material parameters /3,7- Then, in the 
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presence of the dislocation we get: 

e a = % + d a u» + 0{ 7 ), 


= $ - 9^u a 


0 ( 7 ), 


( 52 ) 


a,/x = l,2,3 eft = (1 + 0 ( 7 ))$ 

- 7 . - ' -w/> a,//= 1 , 2 ,3 e“= (1 + 0 ( 7 ))$ 

e = |det e“| = (1 — djU* + 0 ( 7 )), * = 1,2,3 

Af = +V,(u-K W) + 0(J3), A 0 = 0((3) 

This form of the emergent fields allows us to rewrite the partition function as follows: 

Z~ J DTDI'exp^i ^2 ^2 J d3xdte (id^ ± + A± - S M ^ Si± + (h.c.) / (53) 

Let us introduce the left - and the right - handed currents Jj,Jr. They are related to the ordinary 
electric current J and the chiral current J 5 as J = Jl + Jr, J 5 = Jr — Jl- Anomaly results in the non 
- conservation of these currents. Following [4t| we represent anomaly as a result of the spectral flow that 
follows from the explicit solution of Pauli equation. We extend the consideration of j4G] to the case, when 
both emergent vielbein and emergent gauge field are present, and restrict ourselves to the case, when the 
emergent vielbein does not depend on time while e l 0 = 0 for * = 1,2, 3. The expressions for the anomalies in 
the quasiparticle currents have the form: 


(9vJr±) — 

<«±> = 


1 


— d p (B v - At) d p (B a - At), 
-_ 1 _ d ^ Bv _ A ± } dp{Ba _ A ± } 


(54) 


Here + stands for the quasiparticle current at while — stands for that of — K^. This expression does 
not contain the Nieh - Yan term \nvpo g^ e a d p e b a r]ab (that is expressed through the ultraviolet cutoflfl A) 
simply because this term vanishes when the vielbein does not depend on time and e l 0 = 0 for i = 1, 2, 3. We 
do not discuss the possible appearance of this term in the case of the time - depending vielbein. 

Eq. (|54[) m ay be derived in the absence of torsion in a number of different ways (see, for example, 
33, 40, 4ll. l42(| 4. At the dislocation the emergent magnetic field is of the order of 4, where S is the area of 


the dislocation core. Torsion magnetic field is of the order of ^ and is suppressed at low energies compared 
to the emergent magnetic field due to the cutoff A in the denominator. Therefore, torsion may be considered 
as a perturbation. Moreover, in the absence of torsion electric field Tg) it can be shown, that this perturbation 
due to torsion that appears due to the dislocation does not affect the spectral flow and the total expression 
for the chiral anomaly of Eq. (1541) (at least, for the non - marginal cases, when the emergent magnetic field 
does not vanish and therefore dominates over the torsion magnetic field). In Appendix B we illustrate this 
pattern using the direct solutions of Pauli equation in the presence of the torsion magnetic field of particular 
form - when it is directed along the emergent magnetic field. This is the case realized in the presence of the 
dislocation. 

In the marginal case bK = 0 the emergent magnetic field carried by the dislocation vanishes, and torsion 
cannot be considered as perturbation (see Sect. 16.61) of Appendix B). We do not work out this case in details 
and do not exclude, that in this case the emergent relativistic invariance may be broken by the anomaly. 

Since the emergent 17(1) fields at the two Dirac points have opposite charges, and its time derivative 
vanishes, we arrive at the total chiral anomaly (J 5 = Y^±( d R ±~ d L ±)) : 


(«> = 


^- 2 d p B v d p B„ 


(55) 


It is given by the external electromagnetic field B only. The emergent magnetic field contributes the expres¬ 
sion for the anomalous transfer of chirality between the two Dirac points: 


1 


<0„(J£ + - J£_)> = --^r°d p B„d p A t 


(56) 


3 Very roughly it denotes the scale, at which the quasiparticle dispersion loses its Dirac cone form. That means, that A < |K|. 
Here K is the position of the Dirac point, that obeys |K| ~ 1 /a, where a is the interatomic distance. 
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(Here we substitute A = A + ). 

5.2. Anomalous production of quasiparticles at the dislocation 

Now we are in the position to calculate the anomalous production of quasiparticles at the dislocation in 
the presence of external electric field. We suppose for simplicity, that the external magnetic field is absent. 
Therefore, the mentioned above marginal case Kb = 0 is trivial: there is no production of quasiparticles 
at all. At Kb / 0 we consider the gauge, in which B = 0, Bq = <j>. According to Eq. (l54l) the anomalous 
production of quasiparticles appears at each Dirac point given by 

(8^,±) * J o b a K a d^dY%s)S^(X-Y(s)), (d^ ± ) « b a K a d^dY^(s)6^(X-Y(s)) 

(57) 

Here ± refers to the Dirac point. One can see, that the total chirality is conserved = (d^Jf) = 0. 

At the same time there is the transition of chirality between the two Dirac points. 

Let us denote the density along the dislocation of the total production rate of the left/right - handed 
quasiparticles (incident at the two Dirac points iK^ 0 ^) by q R R = 4? w w h ere Ql r * s the total number 

of the left-handed/right-handed quasiparticles with momenta around ±Kl°l. We have: 

^“ ± i( K '‘ bi )( E ‘ 4 *)' » ^ 3/5 ( K , bi ) (E t n“) (58) 

Here the unity vector along the dislocation is denoted by n (we assume n 2 = 1). This expression is valid 
also in the original reference frame (in which the Fermi velocity is anisotropic). 


5.3. Chiral magnetic effect along the dislocation 

The production of quasiparticles of Eq. (1581) in the presence of external electric field results in the 
appearance of the chiral imbalance at each Dirac point. Let us suppose for definiteness, that the dislocation 
is a straight line connecting boundaries while the sample area in the plane orthogonal to the dislocation is 

S±. 

Suppose, that the considered above process of the quasiparticle pumping from vacuum takes place without 
external magnetic field Hr due to the emergent magnetic field Ha incident at the dislocation. The created 
particles and holes leave the dislocation and take part in the scattering processes within the whole volume of 
the semimetal. Following [51| we consider the simplified equation that describes the kinetic of the appearance 
of chiral imbalance at the Dirac point ±K in the presence of external electric field E directed along the 
dislocation: 


pf = ± 


1 


Pt 


A 20 KhE ~ 

47 Ty 


(59) 


Here the chirality changing scattering time is denoted by Ty. while the chiral charge density (averaged over 
the whole volume of the sample) is 2pf = P%~Pr, while pj R is the density of the left-handed/right - handed 
particles at the Dirac point ±K. We neglect scattering processes with the possible transitions between the 
fermions incident at the different Dirac points. The solution of this equation for r>ry gives 


P 5 


± 


4:TT 2 S 


-KhEry 


(60) 


According to 52J in the presence of temperature T and ordinary chemical potential p. the induced chiral 


chemical potential is related to p 5 as follows: 


pf v F = ^2 (^) 3 + + £ 2 ) 


(61) 


(The thermalization occurs mostly out of the dislocation. Therefore, Eq. (43) of (52y can always be applied.) 
For p,T < this results in pf = vp (Sf 2 p^j ■ However, following [5l| we imply the more realistic 
situation, when p, T p 5 , and get 


Pf = 


3 pfvp 
T 2 + K 


(62) 


14 







The quasiparticle currents along the dislocation are given by 


f^(*o, Py , Pz \ee k a a a \* 0 , Py , Pz ) = / o ( K i b ^)dY k (s)6^(X -Y(s)) (63) 

^^(^p v ,pA ee > a \^ Py , Pz ) = ±||r/ 0 (KjyJd^WiPi^-yw) 

We used here, that (e e()n a ) = n^ l (l — du) = n ^ 1 along the dislocation. These currents are directed along the 
dislocation. The net electric current j’cme = Y1±Ur + Jl) originated from the chiral magnetic effect is the 
sum of the currents at the both Dirac points ±K: 



Jcme 


3v|(Kb y E tv 

4Sj_7t 4 (T 2 + £) 


dY k (s)S (3) (X -Y(s)) 


(64) 


The total electric current through the sample is directed along the dislocation and is given by J = aohmE + 
(Jcme-Bj where <rohm is the ordinary Ohmic conductivity while the contribution to conductivity originated 
from the chiral magnetic effect is given by: 


0"CME 



4Sltt 4 (T 2 + £) 


(65) 


6. Conclusions 

Let us summarize the results obtained in the present paper. 

1. Quasiparticles in Dirac semimetal in the presence of dislocations experience emergent magnetic field 
and emergent vielbein. Those emergent fields are expressed linearly through elastic deformations. 
Tensors of material parameters /?, 7 enter the corresponding expressions. 

2. When the material parameters may be neglected, the dependence of emergent magnetic field and 
emergent vielbein on the elastic deformation in the presence of dislocations remains nontrivial. In this 
case the dependence is of pure geometrical origin. It results in the appearance of emergent magnetic 
flux and torsion singularity localized at the dislocation. 

3. The emergent magnetic flux and the torsion singularity carried by the dislocation may be observed 
using scattering of quasiparticles on the dislocation. Emergent magnetic flux results in the Aharonov 
- Bohm effect while the torsion singularity manifests itself in Stodolsky effect. 

4. In the presence of external electric field the quasiparticles are pumped from vacuum. This results in 
the appearance of chiral anomaly at each Dirac point and at the transition of chirality between the 
two Dirac points. 

5. We calculate the anomaly based production of quasiparticles along the dislocation in the presence of 
external electric field. 

6 . The production of quasiparticles along the dislocation results in chiral imbalance. It drives the chiral 
magnetic effect, that is the appearance of electric current along the dislocation without any external 
magnetic field due to the emergent magnetic field incident at the dislocation. 

The given results confirm that the low energy effective theory of quasiparticles in Dirac semimetal (in¬ 
cluding the chiral anomaly) simulates relativistic high energy theory when the Coulomb interaction between 
the quasiparticles is neglected. Although the latter interaction is strong due to the large value of effective 
coupling constant (EsJ, the experience of the investigation of graphene [2] indicates, that for the consid¬ 
eration of various emergent relativistic effects it may be possible to neglect the interactions between the 
quasiparticles. 

We observe, that the electric charge of the quasiparticles in the considered effective theory is conserved, 
while the chirality at each Dirac point is not. Chiral anomaly drives the production of quasiparticles at 
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each Dirac point in the presence of parallel magnetic and electric fields. The dislocations give an effective 
gravitational background for the given low energy theory and the emergent magnetic flux located at the 
dislocation. In this particular case the gravitational background may be considered as perturbation. The 
chiral anomaly is dominated by the contribution of emergent magnetic field and results in the production of 
quasiparticles at the dislocations in the presence of external electric field without external magnetic field. 

Moreover, the analogue of the chiral magnetic effect exists: there is the electric current directed along the 
dislocation (that is along the emergent magnetic field inside the dislocation). This current is proportional 
to the chiral chemical potential that appears due to the mentioned above pumping of the quasiparticles at 
the dislocatioifi. 

The discussed emergent relativistic effects in Dirac semimetals related to emergent gravity and emergent 
gauge field should affect transport phenomena and, therefore, may, in principle, be measured experimentally 
(see also ( 39 I . 53 . 541). Above we already mentioned one of the transport phenomena that is intimately 


connected with the emergent geometry: the appearance of the contribution to conductivity related to the 
chiral magnetic effect, driven by the emergent magnetic field incident at the dislocation. 

The investigation of Dirac semimetal in the presence of the dislocations opens the possibility to observe in 
laboratory the relativistic quantum field theory effects that may take place in particle physics at the energies 
that are not accessible at the present moment. These are the effects of the hypothetical space - time torsion. 
One of such effects (the Stodolsky effect) was discussed in the present paper. 

In the present paper we did not considered the case of the time - depending strain that would result in the 
time depending vielbein. The appearance of the Nieh - Yan term in chiral anomaly in the Dirac semimetal 
(and the related observable effects) in this case remains an open question. It would also be important to 
investigate the influence on the discussed phenomena of Coulomb interaction between the quasiparticles. 
Since this interaction is strong [Hij . most likely, such a study should include numerical lattice methods. 
Finally, we would like to notice, that the constructions discussed in the present paper may be extended to 
the consideration of Weyl semimetals, which, presumably, have a more rich phenomenology. 

The author is greatful to G.E.Volovik and V.A.Miransky for useful discussions. The work is supported 
by the Natural Sciences and Engineering Research Council of Canada and by grant RFBR 14-02-01261. 


Appendix A 


Here we give the solution of elasticity equation Eq. (1531) for the case, when the origin of the dislocation 
is the infinite straight line ([HI, chapter 27, problem 4). We imply the decomposition of Eq. (1511) . 

In the presence of the dislocation we obtain the equation for the continuous part of u l : 


0 = (A <f fc + 


1 


1-2(7 
12 = -i 21 = 1 


V i V fc 


)(u k CODt -b l i kl ^^) +2b l i il 6V\X), 




= i 22 = 0, X 1 = 


X 1 


i*r 


Here we place the origin of the dislocation at X 1 = 0. We come to the following solution: 


t (X) = b l i 


m 

2 R 

2 ir 


l-kl ( lo § 2R + 7 ) 


1 


~A l 2(1 -a) 


1 V*V 


irjk 


-)2 b l i ll 5 (2 \X) 


b l i kl n \X\ 

'^og 2R + ^ + —a 


1 V ' V ^( l 0 g S + ^ 1)|X ' 2 


( 66 ) 


(67) 


The inverse Laplace operator A 1 is to be defined taking into account boundary conditions and the finite 
size of graphene sample. For the sample of linear size R we come to the following expressions: 


A 


-1 


d 2 K e iKaxa 1 ,, \X\ 


J .Y 


(2t r) 2 |A| 2 


-s (loe M +1> ' |X|<<R 


( 68 ) 


4 Ordinary chiral magnetic effect in Dirac semimetals caused by the external rather than emergent magnetic field was discussed 
recently in | 5 l| . 
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(7 is the Euler constant), and 



d 2 K e iRaxa 
( 2^)2 |/\ | 4 


IF 

47r 




|X| <1? 


Up to an irrelevant constant we have 


tPO 


1-2 a 

2(1-^) 


,,. w 1 . |X|e^ 

b l —log- 1 —— 
2tt B 2R 


+ 6 i 


l;il 


1 WX* 

1 — cr 47T 


(69) 


(70) 


Appendix B 

6.1. Pauli equation in the presence of the torsion magnetic field directed along the dislocation. 

In principle, torsion may be considered as a perturbation while the zero approximation corresponds to 
the solution of Pauli equation with flat vielbein in the presence of emergent (or real) magnetic field and 
external electric field. One can check that the perturbation due to torsion does not affect the spectral flow 
at least, when there is the emergent magnetic field and emergent torsion magnetic field that both are caused 
by dislocations (see below, Sect. 16.51) . Below we illustrate this pattern by the solution of Pauli equation in 
the presence of torsion field X)“- directed along the dislocation, i.e. when the dislocation is directed along the 
third axis, the only nonzero component of torsion is Tf 2 . This is the case, realized inside the dislocation, 
where both emergent magnetic field and emergent torsion are directed along the dislocation. For the purpose 
of illustration we shall consider the further simplification with the only nonzero component of torsion T 3 2 . 
The extension to the case, when the components Tf 2 and/or Tf 2 are nonzero is straightforward. 

Thus we consider the case when constant emergent magnetic field Hj and real magnetic fields Hb : as 
well as the electric field E = Eb are directed along the third axis. The external field is considered in the 
gauge B 0 = — zE , while the vielbein has the form: 

e“ = s; + S^S^ThX 1 (71) 

As it was mentioned above, the nonvanishing component of torsion is given by T 3 2 = Th- It corresponds 
to the so - called torsion magnetic field 0- Pauli equation for the left - handed/right - handed fermions 
receives the form (we denote H = + Hb) (40| : 

idt^L.R = ri^’L.R (72) 

The corresponding Hamiltonian has the form 

U^r = -Ez T (a 1 (-id!) + a 2 (p 2 - T H X\-id 3 ) + HX 1 ) + a 3 (-id 3 )) (73) 

It coincides with the usual Hamiltonian for the particle in the presence of electric held E and magnetic held 
H — H — ThP 3 ■ The spectrum of operator Eq. (1731) for E = 0 is given by: 


c(L) 

°l,Sz,P3 

e(L) 

■' 0 ,s 2 ,P 3 


= ± \/P3 + 


H - T H p 3 


{21 + 1) + 2(H-T hP3 )s z , 1 = 1,2,..., a* = ±1/2 


= sign 


(t hP 3 - H ) 


P 3, 


1 

s z = ^sign 


(r H p 3 - H ) 


The wave function of the lowest state has the form: 

ip2V+ip3Z-\ 

^o. Sz , P 2 ,p 3 ,t[x,y,z) = e 


—ThP 3 ~\~H 


_£2_/ 

V ^-T H V 3 + H J 




(74) 


(75) 


Here |s 2 ) is the eigenstate of ^a 3 corresponding to the eigenvalue s z . We omit the normalization factor and 
denote X^ = (X°, A 1 , A 2 , X 3 ) = (t, x, y, z ). In the same way we find the spectrum of the Hamiltonian for 
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the right - handed particles. It differs from that of the left - handed ones by the expression for the lowest 
Landau level (LLL): 


y(R) 

°0,S Z ,P3 




S z 




(76) 


The degeneracy of the level is v(p 3 ) = 


T H p 3 - H 


S _l, where S± is the area of the system in the plane 


orthogonal to the magnetic field. In the following by the subscript || we denote the component of the vector 
parallel to magnetic field. 


6.2. Evolution in time of the wave function 

Evolution in time of the wave function is given by: 

'f> L ’ R (t, x, y, z) = t T L ’ fl (0, x, y, z ) (77) 

In general case this evolution includes transitions between the eigenstates of However, when the values 

of P 3 and t are sufficiently small and we may neglect Th{p 3 + Et ) compared to 77, then the time evolution 
of the wave function for the LLL is given by 


4' 


L,R 

0.P2.P3 


0 t,x,y,z ) « e 


iEtz±sign(H-T H p 3 )i f\p3+Et)dt+ip 2 y+ip3Z-Tj\H-p 3 T H \(x+ H _ V rf ) 1 . . .. 

0 v H3) | - -sign(77 - T H P 3 )) 

(78) 

It gives the solution of Pauli equation for sufficiently small p 3 and t. Notice, that the wave function of 
Eq. (l78l) is not the eigenstate of 'M.^r gi yen by Eq. (l73ll . However, it corresponds to the definite values of 
momenta: 


Py^of 2 . P3 (*> x , V, z ) = P^of 2}P3 (t, x , y, z), 


/• vTf L i R 

P z '^0,p2,p 3 


( t,x,y,z ) = {p 3 + Et)^o’ R p 3 (t,x,y, z) 


Moreover, we have 


V-l!r^o'p 2 , P3 (*) V’ z ) = =Fsign(7J - T H p 3 ) (p 3 + Et)^ R p3 {t, x, y, z) 
where T-L^r is the Hamiltonian for E = 0 given by 

^ { l,r = T (^(-idi) + a 2 ((-id 2 ) - T H X\-id 3 ) + HX 1 ) + a 3 {-ids )) 


(79) 


(80) 


(81) 


The mentioned above extension to the case of nonzero If 2 results in the substitution T^P 3 —> X ) a =2 3 ^hPo. ■ 
We do not need to consider the case of nonzero T-/ 2 because the x axis may always be chosen directed in 
such a way, that Tf 2 = 0. One can see, that inclusion of nonzero Tf 2 does not affect the main property of the 
spectrum: the state of Eq. (1751) with the given value of P 3 carries the definite time - depending momentum 
Pz{t) = P3 + Et and the definite eigenvalue of the unperturbed Hamiltonian £{t) = q=sign (H)p z {t) (for 
sufficiently small p 3 and t). For the values of momenta of the order of H/Tjj , however, the evolution of the 
wave function in time is more complicated. It involves the transitions between different eigenstates of 77*- 0 -*. 


6.3. LLL in the case 0 

We consider the case, when the torsion electric field Tfo vanishes while the torsion magnetic field Th 
is nonzero. We turn on electric field E = Eb and magnetic fields H = Hb — Ha- In the case of interest 
related to the emergent magnetic field located at the dislocation the typical values of magnetic field are of 
the order of 27r/5'disiocationj where ^dislocation is the area of the dislocation core (except for the case, when 
Kb = 0 so that the emergent magnetic field is absent at all). At the same time Th ~ a/Tjdisiocation, where 
a is the interatomic distance. Therefore, H/Th ~ 2ir /a ~ A, and the left - handed states correspond to 
the left - moving particles while the right - handed states correspond to the right - moving particles (or 
vice versa depending on the sign of 77). The given pattern is illustrated by Fig. [I] The situation, when 
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Figure 1: The flow of the eigenvalues of H near the Dirac point in the presence of external electric field E, magnetic field 
H and the torsion - magnetic field T = Th > 0. We represent schematically the dependence of the Hamiltonian eigenvalue £ 
carried by the solution Eq. £78} on ps and on time. (Here ps is the difference between the third component of momentum and 
the third component of the Dirac point). Crossing of the left - handed and the right - handed branches occurs at ps = —Et. 
At ps = H/Th — Et there is the discontinuity of the spectrum for each chirality. In vacuum the states with £ < 0 are 
occupied. Right - handed states are denoted by solid and dotted lines while the left - handed states are denoted by dashed and 
dashed - dotted lines. The value of ps belongs to the interval [—A, +A], where A ~ 2iy/ a. Since Ha ~ 27r/5di s iocation while 
T h cl/ ^dislocation ? the discontinuities of the spectra do not appear within the Brillouin zone if Ha 7 ^ 0, and may be neglected. 


electric field is directed along the third axis of the Brillouin zone (then the coordinates of the Dirac points 
are (0, 0, ±|K^ 0 ^|)) is also illustrated by Fig. [2] 

Both LLL branches of are parametrized by two parameters p3 and p2 . They have the meaning of 
momenta along the 2 axis and y axis only near the Dirac points, where, in addition, the dispersion S(p 2 ,P 3 ) 
does not depend on p2 . We can always organize the parametrization in such a way, that the line p2 = 0 
connects the two Dirac points ±Kl°h For each value of P3 we denote the number of states that belong to 
the given branch and are parametrized by P 2 by v^ R ’ L \p 3 ). If the emergent magnetic field dominates, the 
LLL branch that represents the left - handed Weyl fermion located at K 1 ^ represents the left - handed Weyl 
fermion at — K( 0 ). The total number of states on the given branch for the values of P 3 varying between the 
positions of the zeros of £ is equal to 

N -= f ^^v {R) (P 3 )o(£ {R Ho,P 3 )), N + = j ^Elv( L )(p 3 )e(£( L )( 0 ,p 3 j) (82) 

By L || we denote the size of the system in the direction of magnetic field. Here £^ R \p 2 ,P 3 ) and £^{p 2 ,P 3 ) 
are the dispersions of the two branches of LLL that correspond to the right - handed and the left - handed 
particles at K correspondingly. The band structure in semimetal and our parametrization are organized 
in such a way, that the sign of £ ( - R \p2,P3 ) (and of £^ L \p2,P3)) does not depend on p 2 . iV_ corresponds 
to the shortest path within the Brillouin zone, while N + corresponds to the path that goes through the 
boundary. We denote by v^ (^ 3 ) the number of states with the given value of P 3 of the branch of spectrum 
that corresponds to the right - handed particles at K. Correspondingly iW') (p 3 ) is the number of states of 
the branch of spectrum that corresponds to the left - handed particles at K. 

6.f. Dynamics of the spectral flow 

There are two distinct processes in the presence of electric field. The evolution in time of the states 
with the values of P 3 near the Dirac points is given by Eq. dj. This corresponds to the appearance of the 
quasiparticles/lioles due to the spectral flow. The second process corresponds to the evolution in time of 
the states with the values of P3 away from the Dirac points. Those states may transform to each other as 
it was mentioned above after Eq. m if Th ^ 0. Besides, the transitions may occur also within the other 
branches of spectrum and between the different branches. However, the latter transitions are suppressed by 
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Figure 2: Schematic pattern of LLL in the presence of external electric field E, magnetic field H and the torsion - magnetic field 
Tjf > 0. We represent the eigenvalue 8 of 7^°) as a function of £> 3 , where p 3 is the third component of the total momentum 
(rather than the difference between the momentum and the position of the Dirac point). It is supposed, that the third axis is 
directed along K. Crossing of the branches occurs at the two Dirac points d=|K|. When the emergent magnetic field dominates, 
the left - handed branch at K is at the same time the left - handed branch at —K. (The lines are closed through the boundary 
of the Brillouin zone.) The values of 8 carried by the occupied states at t ^ 0 are shifted by Et (solid lines). As a result for 
E ^ 0 we observe the production of quasiparticles at each Dirac point. 


the difference in energy. Thus, we assume, that the states of LLL with p 3 between the positions of — K and 
K remain occupied. 

The state that corresponds to a certain value of p 3 near the Dirac point changes its energy as £(t,p 3 ) = 
£(0,p 3 + Et). As a result the system leaves its vacuum state (in vacuum all states below £ = 0 are occupied). 
The production of the quasiparticles is given by 


pEt 

(dj R ) = d t 

Jo 


dp 3 v^{p 3 ) Ev( r \ 0 ) 


2ir S± 


2 irS± 


47T 2 


HE, 


» -j^ HE 


(83) 


We see, that the presence of the torsion magnetic field Tf 2 does not affect the expression for the anomaly. 
Our consideration may easily be extended to the case of nonvanishing components T 32 and Tf 2 , that also do 
not contibute to anomaly. 

We come to (dJ±) = 0 and 


{dJ 6 ,±) 


2n 2 


(H b - Hf)E 


(84) 


Eq. (18311 may be easily generalized to the case of the slow varying fields with arbitrary components thus giving 
rise to Eq. (1M1) . We expect that this anomalous production of the quasiparticles should affect transport 
properties of semimetal (see also 39|, 53, 54|). 


6.5. Perturbations due to the arbitrary components of the vielbein 

The extension to the case of nonzero T“ 3 and T 23 and nonzero (with i,j = 1, 2,3) is more involved. 

The corresponding components of the vielbein may be considered as perturbations 8e l a = e* — S l a . The 
correction to the energy of the given state (for sufficiently small p 3 ) gives 

£l,r ~ Tsign(if) (d)p (85) 

We denote d = (0, <5e§, 1 — <5e} —Se 2 ). Taking the leading terms in the expansion in the powers of momenta we 
may consider averages (e“) over the state with p 2 = p 3 = 0 as constants. Boundary between the states with 
positive energies and states with negative energies is the line that contains the point p = 0 and is orthogonal 
to (d). If we neglect the transitions between the eigenstates of with different values of P 2 ,P 3 , then the 
spectral flow results in the shift of this boundary line towards the direction of the z - axis by Et. This gives 
the same production rate of quasiparticles as in the case of vanishing 8ef. 
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Notice, that the correction to the hamiltonian also results in the appearance of the transitions between 
the states with different values of P 2 ,P 3 - Matrix element between the states with different values of P 2,3 is 

{p'lUlp) ~ T^sign(i?)(p + p , )(p'|d|p) =F |(l Vi P) (86) 

As before, we neglect transitions between different Landau levels. In principle, the transitions due to Eq. ABED 
complicate the consideration. However, both mentioned correction to energy and the transitions between 
the different eigenstates of the unperturbed hamiltonian are absent for the case of the emergent fields caused 
by the dislocation (because in this case 5e% = 0). 

6.6. Marginal case H/Th <C A 

As it was mentioned above, there exists the marginal case, when bK = 0 so that the emergent magnetic 
field carried by the dislocation vanishes. At the same time Tfj ~ b a , and we may have H/Th <C A, so that the 
pattern discussed above is changed. (Say, for H = 0 all states in vacuum are left - handed or right - handed 
depending on the sign of Tjj )■ As a result, the production of the quasiparticles with low energies near the 
Dirac points given by Eq. (1551) should be followed by the additional appearance of the electrons/holes with 
energies ~ H/Th- The corresponding states would experience the jump of energy ~ H/Th- This indicates, 
that the given process likely cannot be considered as adiabatic, and the whole consideration fails. For this 
reason, we do not consider the given marginal case in the present paper. 
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